Trends in the elastic response of binary early transition metal nitrides 
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Motivated by an increasing demand for coherent data that can be used for selecting materials with 
properties tailored for specific application requirements, we studied elastic response of nine binary 
early transition metal nitrides (ScN, TiN, VN, YN, ZrN, NbN, LaN, HfN, and TaN) and AIN. In 
particular, single crystal elastic constants. Young's modulus in different crystallographic directions, 
polycrystalline values of shear and Young's moduli, and the elastic anisotropy factor were calculated. 
Additionally, we provide estimates of the third order elastic constants for the ten binary nitrides. 



I. INTRODUCTION 

Nitride compounds are a prominent class of materials 
with applications spanning from protective hard coatings 
(mostly transition metal nitrides, TMNs, of the IIIB-VIB 
group but also e.g. BN or SiN)l'^, to optoelectronic de- 
vices (mostly IIIA and VA groups, but also e.g. ScN or 
TiN)^, to potential hydrogen storage materials such as 

LigNi. 

When it comes to their superior mechanical proper- 
ties such as high hardness and Young's modulus, TMNs, 
and in particular TiN are often the industrial material of 
choice for surface coating of e.g. cutting tools. The prop- 
erties of simple binary compounds can be successfully en- 
hanced by forming metastable alloys, e.g. Tii_a;Ala;Nl, 
which at higher temperatures age-hardens before decom- 
posing into its stable constituents^. Recent studies have 
shown beneficial effects (such as increased oxidation re- 
sistance or retardation of the final decomposition step 
to higher temperatures) of additional alloying elements 
in Tii_j;Ala;N^"^ and other systems^. Another approach 
how to improve material properties is via multilayer de- 
sign, where individual layers are typically simple binary 
or ternary systems^^"^^. 

A modern way designing new and improve current ma- 
terials is to combine experiment with modelling. For sim- 
ple and/or small systems, quantum mechanical first prin- 
ciple approaches can be used. However, these are typi- 
cally limited to several hundreds of atoms and, e.g. multi- 
layers, crack propagation or nanoindentor tip-layer con- 
tact become difficult topics to handle. Here, the contin- 
uum mechanics employing finite element method (FEM) 
proves to be a successful tooli^"— . A key prerequisite to 
perform FEM calculations is the knowledge of the elastic 
properties of the studied materials, which are not always 
experimentally available (e.g., because some phases are 
stable only in the multilayer arrangement, but not as a 
bulk material) . In such cases the elastic constants can be 
provided by the first principle calculations. 

The literature on first principle calculations of early 
TMNs (group IIIB-VB) is vast. The main focus of those 
papers is the electronic structure and related material 
chemistry problems, while the calculation of the elastic 



properties is often only a minor part of the results. Ad- 
ditionally, a lot of those reports focus only on one (or a 
few) systems. There are some exhaustive reports on the 
chemical trends in the early TMNsi^"— , but apart from 
Refs. [13 and [2ll they do not discuss the elastic proper- 
ties. In addition, there are some discrepancies between 
reported values (e.g., Cn for ZrN between 304 GPa^^ and 
616GPa^°) which are worth cross-checking. 

The aim of the present paper is to give a comprehensive 
overview on elastic properties of early transition metal 
nitrides and AIN, since these materials are or have the 
potential to be used as protective coatings^"^. In partic- 
ular, we investigate ten binary systems, AIN, ScN, TiN, 
VN, YN, ZrN, NbN, LaN, HfN, and TaN. We focus on 
the cubic variant (Bl, Fm3m, NaCl prototype), which is 
the stable configuration of all of them apart from AIN, 
NbN, and TaN being metastable in this configuration. 
The single crystal elastic constants are validated by sev- 
eral independent approaches as well as by a comparison 
with available theoretical and experimental data. Subse- 
quently, we calculate directionally resolved Young's mod- 
ulus, anisotropy factors and polycrystalline elastic prop- 
erties of these compounds, and rationalise the trends in 
terms of their electronic structure and bonding. 



II. METHODOLOGY 
A. Deformation modes 

The linear elastic response of cubic materials is fully 
described by three independent components Cxxxx, Cxxyy, 
and Cxyxy of the fourth rank tensor of the second order 
elastic constants (SOECs). It is convenient to represent 
this tensor with a 6 x 6 matrix 
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where Cn — Cxxxx, C12 — Cxxyy, and C44 — Cee — Cxyxy 
Here we make use of the Voigt notation xx ^ 1, yy ^ 2, 
zz 3, yz 4, xz 5, and xy 6. An additional 
relationship links SOECs with the bulk modulus, B, 
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B describes the elastic response of materials to volume 
changes, and it is obtained as a fitting parameter from 
the Birch-Murnaghan equation of state^i. Consequently, 
two other deformation modes are needed to obtain all 
independent components of the cubic elastic tensor. 

The first pair consists of orthorhombic and monoclinic 
deformations. The orthorhombic mode results in a strain 
tensor 
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,0 



(3) 



and the corresponding strain energy density, U{d), is 

UoMS) = £tot{S) - £oq = {Cn - Ci2)S^ + 0{6^) . (4) 

Here, ftot(^) and £cq are the total energies per unit vol- 
ume, corresponding to eorth(<5) and eorth(O), respectively. 
A monoclinic deformation yielding a strain tensor 
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is used to evaluate the C44 elastic constant from the cor- 
responding strain energy density 
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One should note that these two modes keep the unit cell 
volume constant. 

The second pair of deformations, which is also often 
used, is a pair of tetragonal and triclinic distortion. The 
tetragonal deformation corresponds to a strain matrix 
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producing a strain energy density 



(8) 



The C44 elastic constants is obtained from a trigonal dis- 
tortion with a strain tensor 




(9) 



and a strain energy density 

C/tri(<5) = £tot{S) - £oq = 2C44-5' + 0{d^) . (10) 
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5C11 
Cii -f- C12 
§Cii + 3Ci2 

2C*44 

2C44 



6C44 



|Ciii 
iCiii -hC + 112 
§Ciii -I- 3Cii2 -I- C123 
^Ciii -|- 2C144 
^Ciii -I- 2C166 
6C456 



TABLE I. CoefScients from Eq. \T2\ for various deformation 
matrices A1-A6. 



These two deformations are volume non-conserving. 

Recently, Zhao et alM^ proposed a set of six deforma- 
tion matrices which allow for estimation of second and 
third order— elastic constants of cubic materials at the 
same time. These are: 
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The strain energy density in these cases is 

Ua{S) = £tot(<5) - £eq = AS^ + B5^ + 0{5^) (12) 

where the coefhcients A and B are specific combinations 
of Cij and Cijk as given in Table ID 

B. Calculation details 

Quantum mechanical calculations employing density 
functional theory (DFT)21i2^ were carried out using 
Vienna Ab initio Simulation Packag o^^'^" . Projec- 
tor augmented-wave pseudopotentials^ together with 
the generalised gradient approximation (GGA), as 
parametrised by Wang and Perdew'^^, for the exchange 
and correlation potential are used. The plane-wave cut- 
off energies and the fc- vector samplings of the Brillouin 
zone were carefully checked to provide a total energy ac- 
curacy in the order of 1 meV/at or better. They are listed 
in Table HIl together with the used pseudopotentials; the 
suffices _sv and _pv refer to the exact valence configu- 
ration taking into account explicitely also the s and p 
closed-shell electrons, respectivelly. 
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TABLE IL An overview PAW-GGA pseudopotentials, plane 
wave cut-off energies, and the Monkhorst-Pack sampling of 
the Brillouin zone used in this study. 



III. RESULTS AND DISCUSSION 
A. Equilibrium properties 

The optimised lattice constants, a, formation energies, 
Ef, and mass densities, p, are summarised in Table lllll 
Since there have been a vast number of publications on 
experimental as well as calculated equilibrium structure 
parameters of these early TMN compounds (see e.g., 
Refs. d [H, [H, [m, m and S, and references therein), 
we limit the comparison of the here calculated lattice pa- 
rameters to the experimental data from Ref. iBSl . The 
calculated lattice constants are, as expected for GGA, 
slightly larger than the experimental values. The error 
is smaller then 1% except for the case of TaN, where a 
deviation of about 1.5% is obtained. This is likely to 
be related to the fact that cubic TaN is metastable and 
prefers N deficient configurations resulting in a signifi- 
cant decrease of the lattice parameter with respect to a 
stoichiometric configuration^. 

The trends in the energy of formation, i.e. less nega- 
tive values as one moves from the IIIB to the VB group, 
as well as the absolute numbers agree well with those 
presented by Rovere et alM. 

Lastly, from the calculated lattice parameters (equi- 
librium volume) and the atomic weights we computed 
mass densities. Again, apart from the TaN case, where 
the under-stoichiometry of the experimental compound 
is likely to play a role, the agreement with experimental 
values is satisfactory. 



B. Single crystal elastic constants 

When calculating the single crystal elastic constants 
as described in Sec. Ill A[ one should check how the fitted 
CijS depend on the maximum deformation, Smax, i-e., on 
the range of deformations applied to the unit cell. When 
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5.306 
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TABLE III. Calculated lattice constants, o, formation energy, 
Ef, and mass density, p. The experimental lattice constants 
icxp are taken from Ref. [3^, the experimental values of density 
Pcxp are from Ref. [sj and those marked with an asterisks are 
from Ref. [H. 



<Jmax is too Small, the accuracy of CijS is likely to be influ- 
enced by the numerical inaccuracies of the DFT calcula- 
tions, while non-linear elastic (and perhaps also plastic) 
effects are no longer negligible for too large Jmaxr^- To 
illustrate this behaviour, we plot in Fig. [l]the Cn and 
C44 elastic constants of ZrN as a function of (5max and 
the order of the fitting polynomial. It follows, that with 
increasing order of the fitting polynomial, the plateau 
region where the specific elastic constant is independent 
of (Smax, enlarges. At the same time, the onset of the 
plateau shifts to higher values (^max- The reason is that 
a high fitting order leads to an over-fitting of the too few 
data points for a small (Smax- In extreme cases, such over- 
fitting may lead to incorrect plateaus, as shown e.g. for 
C44 (monoclinic deformation) using a 13*^ order fitting 
polynomial. 

In general, the combination of tetragonal and trigo- 
nal deformations gives more robust results for the early 
TMN in the cubic structure then the combination of or- 
thorhombic and monoclinic deformations. Nevertheless, 
the most robust results in terms of plateau values scat- 
ter and the dependence on the order of the fitting poly- 
nomial, were obtained when employing the deformation 
matrices Al and A6 (see Fig.[T]). 

Recently, Udyansky et al.^^ showed that the elastic 
constants of a-Fe are also highly sensitive to the value 
of the smearing parameter, cr, in the Methfessel-Paxton 
scheme. We have therefore checked the convergence of 
CijS also with respect to a. It turns out that the cu- 
bic early TMNs are not hugely sensitive to ti, but in 
some cases, e.g., ZrN or NbN, the elastic constant val- 
ues change by up to 5% when sigma is increased from 
0.02 eV to 0.8 eV. Nevertheless, these variations typically 
take place only for small values of a, and a converged be- 
haviour is obtained for cr p» 0.6 eV. 

The single crystal elastic constants are summarised 
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FIG. 1. (Colour on-line) Dependence of the ZrN elastic constants on the range, — (5max < <5 < (5max, of deformation taken into 
account, and the order of the fitting polynomial. The upper and lower row correspond to Cn and C44, respectively, as obtained 
from various approaches described in Sec. Ill Al 
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FIG. 2. (Colour on-line) Single crystal elastic constants, Cn, 
C12, and C44, and resulting Cauchy pressure, C12 — C44, of 
the compounds investigated here. 



in Table IIVI When possible to evaluate, we give 
the CijS based on all three methods described here 
(i.e., orthorhombic-l-nionoclinic, tetragonal -|-trigonal, 



and deformation modes). Since the yll-l-A6 de- 

formation modes were the only ones to provide well con- 
verged results for all ten binary nitrides, we show them 
in Fig. [2] and we will use them in the following analysis 
for consistency. 

A comparison with previous DFT-GGA literature 
dat a^^i^^i'^^i'^^'^'^i'*'^''^^ yields, apart from a few exceptions, 
a good agreement with our results. The local density ap- 
proximation (LDA) based Cn elastic constants from the 
literature^Sii^'^ are higher than our GGA-based data. 
This is a consequence of over- and under-binding of LDA 
and GGA, respectively, resulting in too small lattice con- 
stants and consequently too hard elastic constants in 
LDA. Finally, although many of the here calculated elas- 
tic constants agree well with the available experimental 
data, in a few cases the discrepancy is as large as 20% 
(e.g. Cii of VN). 

In all cases, the obtained elastic constants fulfil the 
stability criteria for cubic crystals 



C44 > , Cn > \Ci2\ , Cn + 2Ci2 > 



(13) 



The elastic constant Cn is significantly stifFcr than the 
other two elastic constant. Within each row the Cn 
and C12 elastic constants monotonically increase with in- 
creasing atomic number at the same time, Cn decreases 
from Sc to Y to La (isovalent IIIB group) while it in- 
creases from V to Nb to Ta (isovalent VB group). It 
has been suggested in the literatur o^°'^^ , that a negative 
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M, GGA 
13, GGA 
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13, GGA 
3 GGA 
il LDA 
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13, LDA 
23 LDA 
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13, GGA 
S GGA 
43, GGA 



TABLE IV. Single crystal elastic constants, "o" , "t" , and "A" 
stand from values calculated by orthorhombic and monoclinic, 
tetragonal and trigonal, and using the Al and A6 deformation 
modes, respectively, "c" and "e" stand for calculated and ex- 
perimental data from the literature, respectively. The results 
in bold are used for the further analysis of polycrystalline 
elastic properties. 



Cauchy pressure C12 — C44 < corresponds to more di- 
rectional bonding while positive values indicate predom- 
inant metallic bonding. Indeed, the calculated Cauchy 
pressure is most negative for AIN in which a significantly 
larger charge transfer from cation to anion takes place as 
e.g., in TiN'^^. The Cauchy pressure increases to positive 
values with increasing number of valence electrons within 
each periodic table row, as those contribute mainly to the 
metal-metal d-d interactionai^i^. These trends may be 
used in the materials selection process to realize specific 
requirements. 



C. Directional Young's modulus 

The Young's modulus, E, in a certain direction, ^, is 
defined as the ratio of longitudinal stress to longitudinal 
strain in this direction. The elastic compliances, Sij, are 
in the case of cubic crystals the solution of the following 
set of equations^—: 
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where li, I2, and ^3 are the directional cosines of ^ 
(100), (110), and (111) directions this becomes 

-E'/inn\ =1/ Sii , 
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The results are plotted in Fig. [3l The Young's modu- 
lus in (100) follows mostly the same trend as Cn, since 
the Cn elastic constant has the strongest contribution 
to -E{ioo) ■ There is a considerable difference between the 
semiconducting compounds AIN, ScN, YN, and LaN, in 
which the (100) direction becomes the softest, and the 
metallic TiN, VN, ZrN, NbN, HfN, and TaN, where the 
(100) direction is clearly the strongest. In addition, the 
Young's modulus of AIN in the (111) direction is more 
than 1.5-times larger than in any other of the here inves- 
tigated TMNs. This is mainly caused by the high value 
of C44 of AIN, suggesting that AIN is much stronger in 
shear deformation than the other TMN. 

To quantify the anisotropy, we employ the Zener's 
anisotropy ratic^*^, A, defined as 

L-ll — O12 
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FIG. 3. (Colour on-line) Young's modulus in the (100), (110), 
and (111) directions. 
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FIG. 4. (Colour online) Zener's anisotropy ratio, 
£'(111) /^(ioo> J S'lid B/G values for the ten nitrides studied 
here. 



The results, together with the ratio £'(iii)/£'(ioo) which 
provides similar information, are shown in Fig. 21 The 
results suggest that AIN is clearly stiffer in the (111) than 
in the (100) direction. The opposite result is obtained for 
the group IVB and VB TMN where the (100) direction is 
the stiffest. The group IIIB semiconducting TMN exhibit 
values of both, A and i?(iii)/-E(ioo)i very close to 1. This 
implies that their elastic behaviour is almost isotropic. 
The most isotropic response is predicted for YN with 
A = 1.05. The (an)isotropy of the Young's modulus is 
visualised in Fig. [S) Fig. [SJd demonstrates the isotropic 




FIG. 5. (Colour on-line) 3D representation of the directional 
dependence of Young's modulus for (a) AIN, (b) YN, and 
(c) TiN. 



elastic response of YN, while the comparison of Figs. [SJi 
and c shows the qualitatively different elastic response of 
AIN and TiN. 

Some insight into these trends can be gained from con- 
sidering the differences in bonding. The bonds in AIN 
are strongly ionic'^^ while the TMNs contain a significant 
part of the covalent bondingi^. Since the covalent bond is 
stronger than the ionic, this can rationalise why AIN has 
the smallest value of i?(ioo) ■ When going from group IIIB 
to VB elements within each row, the extra electrons fill 
the bonding metal-metal orbitals, while the hybridised 
sp^d? states move to lower energiesi^iS. This can be in- 
terpreted as strengthening of the hybridised sp^d?' bonds 
which are oriented along the (100) directions. As for 
the high C44 elastic constant of AIN, one may argue that 
since there are no d electrons available to form the metal- 
metal bonds in the (110) directions (as it is the case for 
the IVB and VB group elements), upon a shear defor- 
mation a significantly increased repulsion between Al-Al 
and N-N ions as they get closer occurs, which causes the 
high value oi Ei^my 



D. Polycrystalline properties 

Several models exist which assess the isotropic poly- 
crystalline elastic properties using the anisotropic sin- 
gle crystal elastic constants of a given material. Voigt's 
approach^ of constant strains in all grains yields the 
upper limit, Gy and Ev^ to the polycrystalline shear 
and Young's moduli, respectively. On the other hand, 
Reusa^ proposed to apply constant stresses in all grains, 
which yields lower limits and Er. Taking By — 
Bn = B, where the bulk modulus, B, is obtained from 
the Birch-Murnaghan equation of state^^, one gets 



Gv 
Gr = 

Er, — 



Cii — C12 + 3C44 

' 5 ' 
5 

4(5ii - S12) + 3S44 ' 
a — V, R . 



3B + Ga 



(18) 
(19) 

(20) 



Finally, HersheyS^ derived an equation for self- 
consistently calculating the shear modulus, Gh- In this 
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approach, Gh is the real positive root of the following fourth order polynomial 
I 

MG% + 16(4Cii + bCi2)G]j + (3(Cii + 2Ci2)(5Cii + AG12) - S{7Cii - 4Cr2)Cii) Gjj- 

- (29Cii - 20Ci2)(Cii + 2C,2)CmGh - 3(Cii + 2Ci2f{C\i - ^12)^44 = . (21) 

This equation can be simplified by dividing it with (3Cii + 6C12 + 8G) to a third order polynomial-^ with the same 
positive real root 

8G^ + (5Cii + 4Ci2)G2 - G44(7Cii - 4Ci2)G - G44(Cii - Ci2)(Cii + 2C12) = . (22) 



Subsequently, Eq. [20] is used to estimate the Young's 
modulus within Hershey's approach. 

The thus calculated polycrystalline elastic constants 
are shown in Fig.[6]together with the bulk modulus. They 
fit well with the few accessible experimental data-points 
included (black full symbols). The trends in E and G 
are akin: the maximum value in each row of the periodic 
table is obtained for the group IVB TMNs. The spread 
between Gu and Gy (shaded in Fig. HI), as well as be- 
tween Efi and Ey is very small for the group IIIB and 
IVB TMNs, suggesting that the elastic properties of poly- 
crystals of these materials will not be hugely influenced 
by the misorientations of individual grains. A different 
situation is obtained for AIN and group VB TMNs (in 
particular, for NbN and TaN), where the Reuss-Voigt 
range is quite large. As shown in Ref. the ratio be- 
tween the Voigt and Reuss bounds depends non-linearly 
on the anisotropy factor A. The ratio becomes partic- 
ularly large when the anisotropy A approaches 0, as in 
the case of NbN and TaN. As a consequence, these ma- 
terials are expected to be strongly affected by the actual 
microstructures (i.e., not only by the grain orientations, 
but also by the shape of the grain). 

Based on an evaluation of a large experimental data 
set, Pugh^^ proposed that the higher (lower) the B/G 
ratio is, the more ductile (brittle) the material is. This ra- 
tio is plotted in Fig. |4l In general, the ductility increases 
from IIIB to VB group (e.g., with increasing number of 
valence electrons and thus increasing amount of metallic 
bonding), and within each group from lighter to heavier 
elements. 



E. Third order elastic constants 

The methodology employing the deformation matrices 
allows also to easily estimate third order elastic 
constants (TOECs), by following Eq. [T^ and relations 
in Table m TOECs, Cijk, appear in the Taylor series 



expansion of the strain energy 

ij ijk 

= \ Y1 i'^'i + '^vkVk + . . . j mVj , (23) 

ij \ k / 

where £ap = (1 + 5ap)Tji/2 is the relationship between 
components Sap of the Lagrangian strain tensor and 
six components rji of a corresponding vector in Voigt 
notation^. According to the above equation, TOECs 
give corrections when applying such large strains that 
linear elasticity no longer applies. TOECs are thus use- 
ful to describe the pressure dependence of second order 
elastic constants, Cij, or thermal properties of solids^. 
This can be of a particular interest for thin films where 
residual stresses in the range of several GPa can be real- 
ized. 

The computed TOECs for the binary systems inves- 
tigated in this work are summarized in Table [V] For 
isotropic aggregates of cubic crystals, Lubarda^ derived 
equations for Voigt- and Reuss-type averages of TOECs. 
These equations are equivalents to the elastic constants 
expressed by Eqs. [TBVBni The corresponding formulae 
are briefly summarised in Appendix [X] The three poly- 
crystalline Voigt- and Reuss-type TOECs, (7123, ^144, 
and C456 are presented in Fig. [7] as the upper and lower 
boundaries of the shaded areas. These boundaries pro- 
vide an estimate for the expected spread of the data and 
depends on the actual microstructure. The lines in Fig. [7] 
represent the Hill's average 

^ Cijky+Cijk,R , . 

GijkM = ^ (24) 

of the Voigt-type, Cijk,v, and Reuss-type, Cijk,R TOECs. 
The six TOECs describing a crystal with the cubic sym- 
metry are in the isotropic case related by 

Cm = 6*123 + 6G144 -f- 8(7456 , (25) 

^112 — G123 + ^244 , (26) 

Ciee = ^144 -|- 2C456 ■ (27) 
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FIG. 6. (Colour on-line) Calculated bulk modulus, B, and 
polycrystalline Young's, E, and shear, G, moduli. The shaded 
area correspond to the Reuss's and Voigt's limit cases. The 
full symbols denote experimental values from literature. 



The results suggest that in each row of the periodic 
table, Ci23 and C456 decrease to more negative values 
with increasing number of valence electrons from 3 to 5. 
For C144 no clear trend is observed. Since the TOECs 
are mostly negative, second order elastic constants get 
stiffer with compressive stresses while they soften under 
tension (compare with Eq. 



IV. CONCLUSIONS 

Calculating elastic properties using density functional 
theory is a powerful technique, in particular when ma- 
terial phases single crystals are not experimentally ac- 
cessible. In this paper we provided a coherent descrip- 
tion of the elastic behaviour of nine binary early transi- 
tion metal nitrides (ScN, TiN, VN, YN, ZrN, NbN, LaN, 
HfN, and TaN) and AIN. These binary compounds are 
of high technological interest for designing materials with 





Cm 


C112 


C123 


C144 


C166 


C456 




[GPa] 


[GPa] 


[GPa] 


[GPa] 


[GPa] 


[GPa] 


AIN 


-5200 


-400 


330 


320 


-850 


380 


ScN 


-5100 


-190 


260 


200 


-330 


285 


TiN 


-7100 


-370 


430 


175 


-475 


80 


VN 


-8000 


-400 


420 


541 


-450 


-300 


YN 


-4100 


-160 


180 


180 


-225 


250 


ZrN 


-6450 


-310 


370 


150 


-370 


-5 


NbN 


-8600 


-190 


115 


300 


-480 


-235 


LaN 


-1200 


-550 


650 


140 


-80 


180 


HfN 


-7050 


-350 


520 


170 


-450 


-300 


TaN 


-9800 


-20 


-190 


340 


-600 


-310 



TABLE V. Third order elastic constants as obtained from the 
six deformation modes A1-j46 (Eg. I12|) . 




FIG. 7. (Colour on-line) The polycrystalline third order elas- 
tic constants. The shaded are corresponds, for each com- 
pound, to the spread between Voight- and Reuss-type ap- 
proach, while the solid line with data points represents Hill's 
average (Eg. [24)) . 
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application-tailored properties. Single crystal elastic con- 
stants, Cij, and directionally resolved Young's moduli, E, 
in (100), (110), and (111) directions are provided. The 
results clearly indicate the special position of AIN. This 
material has the largest Young's modulus E along (111), 
while all group IVB and VB nitrides exhibit the largest 
E value along (100). These trends could be rationalised 
by analysing the bonding characteristics of these com- 
pounds. Computing the elastic anisotropy we find that 
YN followed by ScN and LaN are the materials closest to 
the elastically isotropic behaviour. Finally, the polycrys- 
talline elastic properties (Young's and shear modulus) 



where A is the anisotropy ratio given by Eq. [TTl It can 
be seen that in case of an isotropic materials {A — 1), 
the two approaches give the same results (compare e.g., 
Eqs.|MlandE2 



were calculated. Good agreement with the rather scarce 
available experimental data was obtained. 

Appendix A: Equations for the polycrystalline 
TOECs 

A detailed derivation of the equations below was given 
by Lubarda— . Here we only summarise the final results 
for the reader's perusal. 

The Voigt type TOECs, Cijky, are obtained from the 
single crystal elastic constants Cijk as 

(Al) 
(A2) 
(A3) 

(A4) 

(A5) 
(A6) 

I 
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C 



123, V 



1 

35 



(Cm + I8C112 + I6C123 - 3OC144 - I2C166 + I6C456) , 



C144 V — 777(^111 + 4Cii2 — 5Ci23 + 19Ci44 + 2C166 ^ I2C456) , 

35 



C456.- 



35 



(Cm — 3Cii2 -I- 2C123 — 9Ci44 -f- 9Ci66 + 9C456) . 



The Reuss-type estimates of TEOCs, Cj* j, ^ , can be calculated as 



1 



456, R — 



hA 



35 V 2A 3 



9 9 
Cm — 3Cii2 + 2C123 ~ ^(^144 — Ciee) + ^C456 



C 



1 



144, R 



A 



2A + 3 



Cm — C123 + -^(Ci44 + 2C166) 



4C 



456, R 



9Ci23,R + I8C1443 + 8C456,R = 9Ci23,V + 18Ci44^v + 8C456,V 
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